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Abstract 

o 

H The Aratyn-Ferreira-Zimerman (AFZ) model is a conformal field theory in three- 

dimensional space. It has solutions that are topological solitons classified by an 
integer-valued Hopf index. There exist infinitely many axial solutions which have 
been found analytically. Axial, knot and linked solitons are found numerically to 
be static solutions using a modified volume preserving flow for Hopf index one to 
eight, allowing for comparison with other Hopf soliton models. Solutions include a 

t-h static trefoil knot at Hopf index five. A one-parameter family of conformal Skyrme- 

Faddeev (CSF) models, consisting of linear combinations of the Nicole and AFZ 

r£3 models, are also investigated numerically. The transition of solutions for Hopf index 

four is mapped across these models. A topological change between linked and axial 
^ H solutions occurs, with fewer models permitting axial solitons than linked solitons 

4^ at Hopf index four. 

1 Introduction 

m 
oo 

The Skyrme-Faddeev model [1][2] is the most well known model which has Hopf soli- 
tions. The model is of physical interest in both quantum field theory and condensed 
matter physics. It has been used as a possible description of potential knotted solitons 
O in multicomponent superconductors [3] [4] and it is of potential relevance in a low energy 

effective theory of QCD [5] [6] [7]. It has a Lagrangian density containing two parts, a 
sigma-type term and a Skyrme-type term, where <f> = (<f>i,<f> 2 , <t>z) maps to the two-sphere: 

S3 

"3 C FN = (d»<f> ■ d^) - K (d»(f> X d"<t>) ■ (d^ X d v <t>) . (1) 

In 3 + 1 dimensions, neither term on its own will allow for the possibility of solitons by 
Derricks theorem [8], but the combination of the terms does. Taking either term to a 
specific fractional power may also allow for the possibility of topological solitons. Taking 
the sigma-term to the power three half's is referred to as the Nicole model [9], 

C Nl = (d»cj>-d,cj>)i . (2) 

The simplest Hopf map is a solution to the Nicole model, and is the only analytic solution 
known. Further axial configurations for higher Hopf charges have been found by imposing 
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axial symmetry [10]. Non-axial solitons in this models have already been investigated 
[11]. Taking the Skyrme-term to the power three quarters is referred to as the AFZ 
model [12] [13], after Aratyn, Ferreira and Zimerman who found infinitely many analytic 
solutions using toroidal coordinates. The Lagrangian density is given by 

C-afz = {Hl„) 1 , where H^ u = cj>- (drf x d u <f>) . (3) 

In these models the target space is an S 2 , requiring that 0-0 = 1. The fractional powers 
for the Nicole and AFZ models are chosen to give scale invariant field theories. 

Each of the above describe a map : IR 3 — >■ 5 2 , with M 3 being topologically equivalent 
to (S 3 , if the field is single- valued at spacial infinity. This is also required for finite energy 
in the Skyrme-Faddeev and Nicole models, and is taken to be 4>{oo) = (0, 0, 1) = e 3 . 
Solitons in these theories will then be classified by an integer valued topological charge 
Q, given by the class of smooth maps (f) : S 3 — > S 2 , i.e. 7i 3 (S 2 ) = Z. This charge is 
topological, meaning that maps of different charges cannot be continuously deformed 
into each other. This charge is referred to as the Hopf index. Pre-images of points of 
the target S 2 are generally closed curves in IR 3 . The Hopf index counts the number of 
times that the pre-image curves of two separate points are linked together. Continuous 
deformations of the inverse map must leave the linking unchanged, so each Hopf index 
splits the set of all maps, <fi, into distinct sets. This means that any static solutions will 
have a given topological charge, Q, and the static solution with the lowest energy for a 
given Q will be a stable soliton. 

The purpose of this paper will be to find the minimal energy solutions for Hopf charges 
one to eight using numerical simulations of the full nonlinear field theory. Minimisation 
will utilise volume preserving flow [11] [14] to overcome technical difficulties relating to 
scale invariant field theories and simulations will not be restricted by symmetry. Fur- 
ther technical difficulties appear in the AFZ model resulting from an infinite-dimensional 
symmetry group acting on target space [15] [16]. The AFZ model has area preserving 
diffeomorphisms of the target S 2 that are symmetries of the equations of motion. Com- 
bining these with the conformal base space symmetry results in a Noether charge which 
has a divergent flux. Examples of this flux diverging along an infinite axis [15] and along 
a ring [17] have been found. Sectors of the model with different Hopf index are also 
not separated by an infinite energy barrier, where a symmetry transformation relates all 
solutions to the vacuum [18]. Solutions to the AFZ model have been found [13] so this 
problem can be overcome (see [18] for details). Fixing the value of <f>(oo) will do this. 
The Nicole model does not have the same target space symmetries [19], so adding a small 
amount of the Nicole model to the minimisation of the AFZ model will break the target 
space symmetries. This will be the approach used below, where a contribution of around 
1% of the energy from the Nicole model will be sufficient. 

The addition of the Nicole symmetry breaking term to the minimisation motivates 
investigation of a one-parameter family of conformal Skyrme-Faddeev (CSF) models. The 
set consists of all linear combinations of the Nicole and AFZ models where the coefficients 
sum to one. The transition of the energy and topological changes of solutions to these 
models will be simulated numerically using volume preserving flow. The investigation 
will concentrate on a Hopf index which has topologically different solutions in the Nicole 
and AFZ models. Details of these models will follow the discussion of the AFZ model. 
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2 AFZ Model 



The interest here will be finding static solutions of the AFZ model, and so this system is 
best described by its static energy, 
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x, = 4>- (di4> x dj<t>) , 



(4) 



with spacial indices, i,j — 1 . . 3. The normalisation gives the charge Q = 1 soliton an 
energy, E AFZ = 1. The position of the solitons will be defined where the field is antipodal 
to the vacuum, i.e. 0=-e 3 . These will be closed string(s) in M 3 . 

It is possible to find an infinite number of solutions to the AFZ model by exploiting 
its symmetries in toroidal coordinates [13]. This leads to an energy that depends only 
on the windings of the soliton about the two cycles of the torus, m and n. The axial 
solitons are denoted A n>m with Hopf index Q = nm. The solutions are then dependent 
on a profile function, which can be found explicitly for any axial soliton. The energy is 



E, 
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2 2 ^/jnjjmj (\n\ + \m\). 



(5) 



This can be re-expressed using a ratio of these winding numbers to show that the energy 
of an axial soliton, for a given Q, is minimised when the two winding numbers are as 
close as possible to one another. 
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The Skyrme-Faddeev model has an energy lower bound of the form Esf > csfQ*, 
where csf is a known constant [20]. No proof for a lower bound of this type exists for 
either the Nicole or AFZ model. If such a bound exists then in the Nicole model the 
maximum possible value for the constant is c^i = l from the exact Q — l solution, where 
the energy is normalised to 1 [9]. The work in [11] suggests solitons for Q > 2 have an 
energy ~ 10 — 12% above this conjectured bound. In the AFZ model the largest possible 
value for the constant is also cafz — 1 from (6) [13]. This bound is attained for all square 
charges Q = n 2 . Clearly, axial solutions when the Hopf index is a prime number will fall 
well above this bound. In fact, static solutions to the Nicole model for Hopf index five 
and seven have an energy in the AFZ model, E^fz ~ 10 — 15% under the energy of the 
axial solutions (6). Although these are not static solutions for the AFZ model, it does 
suggest static non-axial solitons must exist for some Q. 



2.1 Volume Preserving Flow 

The local energy density £afz is given by 

E AFZ = Hl+\ AFZ {l -0-0), where E% = H. (7) 

The Lagrange multiplier, Xafz, constrains the system to the S 2 . Variation of the energy 
density leads to the static equations of motion for the AFZ model, 

H-l [(didrf ■ dj<f>) A - (didj4> ■ d k cf>) (d 3 cf> ■ d k <t>)\ [d l( f> (d t cf> ■ dtft-ditA] (8) 
+H~* [V 2 0A + di<t> {didrf ■ dj<f>)-didj<j> {d l( j) ■ dj4>)-dj<f> (V 2 • d^)] + \ AFZ 4> = 0, 
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where A= (<9j0 ■ 40) and V 2 = <9j<9j0. Xafz is obtained from (8), using 0-0=1. The 
gradient flow equations, with flow time parameter, t', are 

They give the direction of steepest decent of the energy functional (4), where 

f afz = [v 2 0A + d l( j> (440 ■ djft-didrf (40 ■ 40) -40 (v 2 ■ 40)] (io) 

+H~i [(440 ■ 40) A-(440 ■ 40) (40 • 40)] [<9/0 (40 • ^0)-40A] , 

and • F' = is needed to keep the flow on the unit sphere, so \afz = ~4> • F afz- 

An additional symmetry breaking term can to be added to the minimisation to over- 
come technical issues in the AFZ model. The local energy density becomes 

£ = H* +eAi +A(l-0-0), (11) 

again using a Lagrange multiplier A and e is a constant that weights the impact of the 
Nicole model on the flow. Since energies in the Nicole model and AFZ model are of the 
same order, e = 0.01 would add ~ 1% of Nicole model to the overall flow direction. The 
new direction of flow, F, will be given by F = F afz + cF Ni + A0, where the gradient 
flow of the AFZ model, F 'afz, is the same as (10) and the gradient flow of the Nicole 
model, F^i, is 

F Nl = V 2 0A5 + d t 4> (didj<l> ■ 40) A-5. (12) 

Minimisation using the gradient flow, F, will not find stable solutions. Zero modes 
associated with changes of scale will be broken when the theory is placed on a lattice. 
These modes will cause the solitons to shrink until their scale becomes too small and the 
solutions fall through the lattice. This behaviour is seen in the 0(3) sigma model [21] and 
can overcome using a novel lattice configuration [22], but no such lattice configuration is 
known here. This issue has already been overcome for the Nicole model using a modified 
gradient flow that preserves the size of the soliton [11]. This volume preserving flow 
[14] [11] can be used to constrain the minimisation to a direction which will leave a given 
total integral of a given function, v, unchanged. This flow is given by 

This particular flow will preserve, in general, the size of the soliton. Since the soliton 
position is defined to be where 3 = — 1 and </> 3 — > 1 at spacial infinity, this choice of 
volume counting seems reasonable. This flow will reduce the energy, E, and keep the 
volume, V, constant. A simple adaptation of the proof in [14] will show this, but both 
should be clear from the construction of the volume preserving flow. 

Volume preserving flow will be applied to configurations of Hopf index Q — 1..8 using 
initial conditions generated by rational maps [23], details of which will follow shortly. 
These configurations will exist in a finite region, Q, and the field will be fixed at = 
(0, 0, 1) on dfl. For the purpose of numerical investigation Q will be a cubic lattice, with 
unit spacing, of length 150. Derivatives will be approximated to fourth order accuracy 
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using a finite difference method and the fields evolved using an explicit method with 
timestep At < 0.1. The AFZ energy and volume for configurations will be 

Eafz — — l — [ H*d 3 x, V= I vd 3 x. (14) 
167r 2 24 Jn Jn 

These integrals and inner products will be evaluated as summations over the lattice. The 
timestep will be reduced, where appropriate, to allow for continued evolution for a given 
e which may be necessary if symmetry issues of the AFZ model become problematic. 
This will be done to keep e as small as possible. Simulations will use e = 0.01 in general, 
but e < 0.05 may be used for small periods of evolution to circumvent symmetry issues 
which cannot be overcome efficiently by lowering the timestep. In all cases except for the 
charge 1 soliton e > is needed, although minimising a charge 1 soliton with e > reaches 
almost the same static solution as without the symmetry breaking. 



2.2 Initial Conditions 

Initial conditions will be generated using rational maps [23]. The fields are constructed 
by mapping (x\, x 2 , x 3 ) G M 3 to the unit three-sphere using a degree one spherically 
equivariant map. In terms of complex coordinates on the three-sphere, Z\,Zq, with 
|^i| 2 + |^o| 2 = l, 

(Zi,Z ) = ( (xi + ix 2 ) —^-,cosf + ix 3 ^^- J . (15) 



Here, f(r) is a monotonically decreasing function of the radius, r. It has boundary 
conditions /(0) = 7r and /(oo) = 0. As our solutions will exist in a finite region, Q, f = 
on dfl. Initial conditions are then generated when the stereographic projection of is a 
rational function of Z 1 and Z , such that 

1 + 03 q(z u z y liDJ 

where p and q are polynomials in Z 1 and Z , and W is the Riemann sphere coordinate. 
Axial solitons can be found by taking p(Zi) and q(Z ) such that (16) becomes 

yn 

W==^. (17) 

These are A n , m type fields with a Hopf index Q = nm. 
Knotted fields can also be generated from (16) when 

W=7 ^L_ (18) 
Zl + Z b y ! 

These are torus knots denoted K, a ^ and have Hopf index Q = ab+(3a. Here, a > b and are 
coprime positive integers with a > and /3 > 0. The only knotted solution of interest 
here will be the trefoil knot, /C 3 2 . 

Linked configurations can be generated from rational maps with reducible denomina- 
tors, q. Two-component links are denoted C^ b m . The integers n,m represent the Hopf 
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AFZ Energy 


Axial Solutions 




Numerical Results 


Q 


Type 


E 


E/Qi 




Q 


Type 


e 


E 


E/Qi 


1 




1.000 


1.000 




1 


Ai,i 


0.01 


1.003 


1.003 


2 


A 2 ,i 


1.732 


1.030 




2 


A 2 ,i 


0.01 


1.738 


1.033 


3 


A 3>1 


2.449 


1.075 




3 


A 3 ,i 


0.01 


2.414 


1.059 


4 


A 2 ,2 


2.828 


1.000 




4 


A 2 ,2 


0.01 


2.852 


1.008 


4 


A 4:1 


3.162 


1.118 




5 


K 3>2 


0.01 


3.391 


1.014 


5 


A 5 ,i 


3.873 


1.158 




5 


7-1,1 

^1,2 


0.01 


3.451 


1.032 


6 


A 3 , 2 


3.873 


1.010 




6 


A 3 , 2 


0.01 


3.907 


1.019 


7 


A 7 ,i 


5.292 


1.230 




7 


K 3>2 


0.01 


4.407 


1.024 


8 


A 4 , 2 


4.899 


1.030 




8 


A 4 , 2 


0.01 


4.954 


1.041 



Table 1: (left) List of AFZ axial soliton energies, calculated using (5) and (right) are the energies of the 
static solutions found using the volume preserving gradient flow method. The amount of Nicole model 
used to find the solutions, e is also given. 

charges of the individual components and a, b the number of times each component links 
with the other. They will have a Hopf index Q = a + b+n+m. On such example is 

W = ^ = ^ + ^ (19) 

Zl-Zl 2(Z 1 -Z )^ 2(z 1 + z y { ) 

This is an £^ l n link with Hopf index Q = 2n+2 consisting of two ^4 ni i's linked once each. 
Other links can be produced from similar rational maps. 
See [23] for further details of Hopf soliton rational maps. 



2.3 Solitons in the AFZ model 

The previously defined rational maps create field configurations which have a given Hopf 
index. By minimising the energy of different rational maps, potential local and global 
energy minima may be found for a variety of Hopf charges in the AFZ model. This 
section will classify those results, as well as noting transitions between different types of 
configuration. 

The soliton energies for both the axial solutions and numerical results are given in 
table 1 and are plotted in figure 1. The minimal energy soliton positions and linking 
are shown in figure 2. All local minima, as well as details on initial conditions, will be 
discussed. The soliton plots (figure 2) show a level set of (p 3 = —0.9, giving an idea of 
the position of the soliton and a level set 0= (a//x (2 — /i), 0,/i — 1), with \i = 0.1 which 
shows the linking. This is an arbitrary choice of a second pre-image curve used to see the 
linking number, and so the Hopf index, of a given configuration. 

Hopf index one and two solitons are the same as in the Skyrme-Faddeev and Nicole 
models. They are the axial Ai^i and «4 2 ,i configurations, see figure 2. The Ai t i axial 
soliton has an energy of 1.003 for a volume V— 115000. This matches well with the exact 
axial solution and is < 0.5% above the expected energy. The A 2) \ axial soliton has an 
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Figure 1: Plot of E/Q* against charge from table 1 for both the exact axial solution (5) (red), and 
the numerically calculated energies (blue). 

energy of 1.738 with V = 100000. This is also < 0.5% above the energy of the exact 
axial solution. Hopf index three solitons are also the same as in the Skyrme-Faddeev 
and Nicole models, the ^3,1 twisted axial configuration. It has an energy of 2.414 when 
V = 60000. An axial .4.3,1 can be found as a saddle point of the minimisation if held by 
symmetry, with an energy of 2.476 with a volume V — 60000. The axial solution held by 
symmetry is ~ 1% above the exact axial solution while the twisted .4.3,1 soliton is ~ 1.5% 
below the energy of the exact axial solution. The Q = 3 soliton has an energy which is 

3 

around 6% above the conjectured energy bound E > Q*. 

As expected from (5) the minimal energy solution for Hopf index four in the AFZ 
model is an axial ^2,2, see figure 2. This is the same as in the Skyrme-Faddeev model 
but not as in the Nicole model where the solution is an C\'\ link. The ^2,2 solution has 
an energy of 2.852 for V = 130000, < 1% above the expected energy from (5). Initial C\\ 
configurations evolve towards an ^4 2 ,2 type configuration but do not quite make it the 
whole way. The two strings do not seem to be able to unlink themselves and the result 
is a linked configuration but with the two unlinks lying almost on top of each other. 

3 

At Hopf index five the axial .4.5,1 exact solution (5) is ~ 15% above the E > Qi 
conjectured energy bound. Evaluating the AFZ energy of the static configuration in 
the Nicole model gives a much lower result than that of the exact axial solution at 
<2 = 5. Volume preserving minimisation using the symmetry breaking term results in two 
local minima with energies much lower than the ^4 51 axial solution. The lowest energy 
configuration is a ^3,2 trefoil knot, with an energy of 3.339 for a volume V— 150000. This 
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Figure 2: Hopf index Q = 1 — 8 minimal energy soliton solutions in the AFZ model. The numbers 
indicate the Hopf index Q. 

is ~ 1.5% above the conjectured bound. The isosurface near to the soliton's position forms 
a band-like structure rather than the usual string-like structure seen with trefoil knots in 
the Nicole or Skyrme-Faddeev model, (see [11] and [23]). The other static solution is an 
C\ 2 link, like the minimal energy solutions in the Nicole and Skyrme-Faddeev models. 
It has an energy of 3.451 when V = 190000, which is ~ 2% higher than the knotted 
solution. The trefoil knot appears to be the minimal energy soliton at Q = 5, but the 
small difference in energy and the fact that no configuration transitions are seen makes 
this proposition inconclusive. Axial .4.5,1 initial conditions tend to twist, similarly to the 
.4^1 configuration. An approximately axial A§,i configuration with a volume V = 225000 
has an energy of 3.924 which is ~ 1.5% above the expected value. This configuration 
is not static as it tends to twist itself until a section of the string pinches off and the 
solution collapses. Conducting minimisation on a lattice with axial symmetry may result 
in an axial configuration with an energy much closer to the exact solution. 

Hopf index six has a static solution that differs from the Nicole and Skyrme-Faddeev 
models. As the A^,2 exact axial solution (5) has an energy that is only 1% above the 
conjectured bound it seems likely it will be the minimal energy solution. A static A^^ 
configuration has an energy of 3.907 for a volume V — 160000. This is < 1% above the 
energy of the exact axial solution, as was the A22 solution at Hopf index Q = 4. A range 
of initial conditions can be used, where the resultant minimisation leads to an A32 type 
configuration, or something near to an A% t 2- The same issue with configuration transition 
occurs at Q = Q as did at Q = 4. The £2 2' ^i'i' ^3,2 an d ^4,3 initial configurations all move 
towards an .4.3 2 type configuration, but the position strings seem unable to reconfigure 
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themselves into the correct formation. The result are knots or links where the strings all 
approximately in a plane, almost on top of each other. 

Hopf index seven has only one solution and it is the same type as in the Nicole and 
Skyrme-Faddeev models, a trefoil knot. The /C 3i2 has an energy of 4.418 with V— 180000. 

3 

It is ~ 17% lower in energy than the exact axial solution (5) and is ~ 2.5% above Q*. 
All possible starting configurations lead to a trefoil knot when Q — 7, including £^3> £i% 
/C5 2 and /C4 3 type initial configurations. Axial initial conditions twist once again under 
minimisation. In fact, they twist too quickly for the axial profile function to minimise 
and an axial energy to be found. This is likely to be the case for all prime Q > 7, where 
the minimal energy configurations are likely to be knots and links, as in the Nicole and 
Skyrme-Faddeev models. 

3 

The Hopf index eight axial A^2 exact solution (5) has an energy 3% above Q*. The 
minimal energy solution found by volume preserving minimisation from a range of initial 
conditions is a twisted .4.4,2 type configuration, denoted Aa,2 ■ It has an energy of 4.954 

3 

for V = 130000, around 4% above Q*. That makes it ~ 1% above the exact axial 
energy (5). This would be in line with the energies found when Q — 4, 6, except this 
solution is twisted. The configuration seems to be robust, since several types of initial 
conditions lead to an type configuration, including £33, £?2n ^-5,2 an d -4.4,2 type 
fields. As with the axial solutions for Q = 4 and 6, the minimisation of other knot 
and linked configurations results in configurations close to the Aa,2 configuration, where 
the strings are not able to reconnect in the correct orientation under volume preserving 
minimisation with symmetry breaking. Imposing axial symmetry in the initial conditions 
can lead to an Aa,2 axial configuration, which has an energy of 4.978 when V — 110000. 
This is ~ 1.5% above the exact axial solutions energy. As the energy of the twisted 
Aa,2 configuration is approximately what might be expected for the energy of an axial 
configuration under volume preserving minimisation it is not clear whether this is the 
minimal energy Q = 8 soliton configuration in the AFZ model or not. The twisting 
may be an artefact of the Nicole symmetry breaking term in the modified minimisation. 
The twisting may minimise the Nicole component to such an extent that breaking the 
symmetry kills the axial solution. Investigation using a larger domain, Q, may help shed 
light on this. Perhaps another method of breaking the symmetry could also be used to 
back up these results. In either case, the Hopf index eight soliton differs from the solitons 
in the Nicole and Skyrme-Faddeev models. 

2.4 Increasing e 

The results presented above were found using mainly e = 0.01. On occasion short runs 
using e = 0.02 — 0.05 may be performed to overcome breakdowns in the minimisation 
process resulting from symmetry problems in the AFZ model. Taking e > 0.05 during 
minimisation can have a large influence on the resulting static solutions for some Hopf 
charges. Increasing e for Q— 1, 2, 3, 4, 5 and 7 does not alter the topology of the solitons. 
The only change is in the energy, where the field has been slightly altered to minimise 
the increased proportion of Nicole model energy present in the minimisation process. 

At Hopf index 6, minimising the knotted configurations with a larger e leads to a 
A^3,2 type configuration. For small e, they minimise towards a /C 3i2 where the position is 
squashed into a plane and the string sections are close together. For larger e the strings 
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move apart, but the knot still remains somewhat squashed. The linked configurations 
evolve towards an type configuration when e > 0.05. This might be expected since 
the minimal energy solution in the Nicole model is of this type. The Asp configuration 
can also be altered by increasing e. Taking e between 0.05 and 0.10 can evolve the field 
towards a twisted A^,2 type configuration, but with increased AFZ energy. In all cases, 
the energy is significantly larger than the small e minimisations. 

Similar results occur at Hopf index 8. For e < 10 the evolution tends towards an 
A^i twisted configuration. If e is increased further the configuration splits into an C 2 ,' 2 , 
type link configuration. An £3^ type configuration evolves towards an ^4 4 2 twisted 
configuration for low e, but the two strings do not seem to be able to cross over to the 
configuration using the symmetry breaking minimisation. Without the symmetry 
breaking the minimisation breaks down. Increasing e leads to an £33 type configuration 
like that seen in the Skyrme-Faddeev model, where they lie close to each other almost 
in a plane. This is different to the static £33 configuration in the Nicole model, which 
has two unlinks linked perpendicular to each other. Again these configurations have 
a larger AFZ energy than for small e. This behaviour might account for the broken 
symmetry minimisation resulting in an Aa,2 twisted configuration rather than an axial 
Aa,2 configuration for small e. 



3 Solitons in Conformal Skyrme-Faddeev Models 

The Nicole and AFZ models have been investigated both analytically and numerically. 
The most recent numerical developments have been detailed in the previous two chapters. 
Any linear combination of these two models will also result in a scale invariant theory 
with topological solitons. The set of all linear combinations of the Nicole and AFZ model, 
where the coefficients sum to one, form a one parameter family of conformal field theories 
permitting Hopf solitons. They will be referred to as conformal Skyrme-Faddeev (CSF) 
models. 

Differences in the topology of minimal energy solutions at a given Hopf index in the 
Nicole and AFZ models has already been noted. This is clearest in the case of square 
charges, Q = n 2 . Solitons with Q = n 2 in the AFZ model are axial with equal winding 
about the two torus angles. Solitons with Q = n 2 in the Nicole model are links and knots. 
Investigation of solitons in the conformal Skyrme-Faddeev models and the topological 
transition of solutions across the set could shed some light on the nature of the solutions 
of the Nicole and AFZ models. They may also lead to some insights about the Skyrme- 
Faddeev model and its solitons. 

The one parameter family of conformal Skyrme-Faddeev models, £(#), are 

C{6) = cos 2 (9)C Nl + sm 2 (9)£ AFZ . (20) 

Both Lagrangian components have been previously introduced (2), (3), and 6 e [0, |]. 
The two extremes of the set give each model individually, the Nicole model is recovered 
when = and the AFZ model is recovered when 9— |. The static energy of this set of 
models is 



J 327r 2 V2cos 2 (#) + 16tt 2 23 sin 2 (#) 
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This energy has been normalised so that E = Etl = 1 for Q — l, i.e. the energy of the 
simplest Hopf map is set to one in both the Nicole and AFZ models. 

3.1 Volume Preserving Flow 

The CSF models are clearly scale invariant V0, so energy minimisation can be conducted 
using volume preserving flow. The energy densities are 

E e = cos 2 (9)£ Ni + sm 2 (6)£ AFZ + X e (1 -</>•</>) . (22) 

Each model will require its own Lagrange multiplier, Xg and the energy densities of the 
Nicole and AFZ models are as defined previously, i.e. 

£ Ni = (drf ■ = A§, £ AFZ = x d 3 4>) ■ (d t <f> x dj<f>)\* ee iff. (23) 

The equations of motion for this model are just a linear combination of the individual 
equations of motion, leading to gradient flow given by 

F e = cos 2 (6)F m + sm 2 (9)F AFZ + A e 0. (24) 

F Ni , (12), and F AFZ , (10), have previously been given. The Lagrange multiplier is found 
as the flow must map to solutions on the two-sphere, i.e. <fi ■ Fg = 0, such that 

\ = - cos 2 (fl)0 • F Nl - sin 2 (#)0 • F AFZ . (25) 

Volume preserving flow can be applied to fix a scale for the energy minimisation. 

3.2 Solitons and Transitions 

The transition of soliton energy and topology is investigated for Hopf index Q = 4. Static 
configurations at either end of the spectrum differ topologically. This would also be true 
for all square Hopf charges, Q = n 2 , with n > 1 as well as <2 = 5, 6 and possibly 8 (see 
section 2 and [11]). Hopf index Q = 4 is the lowest charge for which this occurs. The 
main questions are: 

• Do configurations alter slowly and smoothly as 6 changes? 

• Do two local minima exist for some values of 61 

• Do two topologically different solutions have the same energy for some 61 

• Does an intermediate soliton configuration exist in between the two extreme 0's? 

• Does a sudden topological change in solutions occur at a given 61 
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Table 2: Static CSF energy for (left) C\\ (right) A2.2 type initial conditions. C l {^ configurations with 
9 > 0.87 move towards an _4 2 ,2 type configuration. A2.2 energies with 9 < 0.87 are the energies of the 
axial configurations, which are saddle points of the minimisation. A2.2 configurations with 9 <C 0.87 
do not have static solutions as the field will attempt to unwind to a single string and collapse. Axial 
configurations for 9 just under 0.87 can re-link and evolve towards the C\'\ minimal energy configuration. 
This transformation takes a very long time, due to the shallow gradient of the restricted energy functional 
in this region. Continued minimisation of the C\'\ solutions for 9 > 0.87 leads towards near-axial 
configurations, with an energy slightly larger than that of the _4 2 ,2 minimal energy solutions. 
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Fi gure 3: Plots of the CSF energies for both configurations for a range of #'s, details found in table 2. 
(top) Full range of 9 and (bottom) region near transition C\\ O A2.2- Plots are of Ej Qi against sin 2 0, 
the coefficient of the AFZ term. 

3.2.1 Variation of 6 

The Nicole model soliton, at Q = 4, is an C\'{ type linked configuration, with no other sta- 
ble local minima existing [11]. In the AFZ model, it has been seen earlier that the soliton 
configuration representing the global energy minimum is an A22 type field configuration. 

The numerical calculations are conducted using the same details as in section 2. All 
solutions will have a volume ^ = 160000, resulting from initial configurations determined 
by rational maps ^2,2 and C\'\ detailed previously. Minimisation where 9 ~ | will require 
a much reduced timestep to circumvent issues relating to the AFZ model discussed in 
section 2. 

Table 2 shows a breakdown of the energies over a range of 9 for both the link and the 
axial solitons. Figure 3 plots these results for the full range of 9 (top) and in the range 
of topological change (bottom). For Hopf index 4, C\\ type configurations exist and are 
very similar in appearance when 9 < |, see figure 4. Axial ^2,2 solitons exist and are 
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Figure 4: Position (red) and linking (green) plots for C\'\ minima for a range of 9. Above 9 — 0.87 

solutions move towards A22 type configurations. For 9 7r/4, C\'\ initial conditions will lead to a 
similar configuration as with 9 — 0.88 and have a slightly larger energy than an A2.2 minima for the 
same 9. 




Figure 5: Position (red) and linking (green) plots for ^4 2 ,2 minima for a range of 9. Configurations with 
9 R3 0.9 fluctuate between an ^2,2 state like 9 = 1.00 and a slightly linked state like figure 4 (9 = 0.88) 
with very little increase/decrease in energy. .4 2i 2 configurations with 9 < 0.86 tend to try to unwind to 
a highly twisted — > C\'\ and collapse along the way as a section becomes thin. 
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similar in appearance for 9 > |, see figure 5. The spacing between models considered, 69, 
is decreased for 9 > ^ due to issues relating to the AFZ model, as discussed in section 
2. The range 9 E (|, |) is somewhat more complicated and more models are considered 
in this region. Selected results are plotted as isosurfaces near the position and linking in 
figures 4 and 5. 

The two components are self intersecting in the Nicole model and as 9 — y | the self 
intersection is broken. For \ <9 < 0.87 the unlinks move inward, eventually towards an 
v4.2,2 type configuration when 9 > 0.88. Transitions between a configuration like figure 4 
(9 = 0.81) and the minimal configuration for a given 9 in this range occurs over a very 
small energy range. For = 0.80, the axial solution is only 0.2% above the linked solution. 
This reduces quickly as 9 — y 0.87. Axial solutions lie in definite minima for 9 J> |. In the 
range 0.87 > 9 > |, solutions fluctuate between an axial configuration and a configuration 
similar to figure 4 (0 = 0.88) during minimisation. The energy fluctuates up and down by 
a fraction of a percent. This indicates that either there are two local minima separated 
by a very low energy barrier, and the minimisation cannot properly distinguish them in 
this range, or that only one shallow local minima exists. Perturbations about this minima 
would then have a very similar energy if the gradient of the energy functional is small for 
a wider range of field configurations. 

Clearly the topological transition of solutions for Q = 4 occurs for 9 > |. This can be 
easily explained by analysing the simplest Hopf map for the Nicole and AFZ models in 
the context of the CSF models. The Q = l energy in the Nicole model is 327r 2 v / 2 ps 467 
and in the AFZ model it is 167r 2 24 ps 266 when they are unsealed. This means that for 
the rescaling of the CSF models, each term gives an equally weighted contribution when 
9 ps 0.91. Minimisation effects would have an equal contribution from both terms for the 
Hopf index one solution here. This matches with the results for the transition seen. One 
might expect the transition of the Q = 4 soliton to occur slightly later than is found here. 
In the AFZ model <2 = 4 is a square charge and so the solution is expected to obtain the 
energy bound 167r 2 23Q4 whereas the Q > 2 solitons are around 11% above the energy 
bound 327r 2 v / 2<2 3 - Note that these bounds are not proven, they are simply the largest 
possible bound for behaviour like 

£>cQs, (26) 

where the constant c is the energy of the simplest Hopf map in each case. These are the 
largest possible values of c allowed. One might then expect a topological transition to 
occur at around 9 ~ 0.94 rather than the 9 ~ 0.87 seen. The cause of the fluctuations in 
this region is the most likely explanation for this. 

4 Conclusion 

Non-axial solitons were found to be static solutions of the AFZ model using numerical 
simulations utilising a modified volume preserving flow. It has been shown that all 
solutions for Q < 9 are axial, twisted axial or knotted solutions. This does not rule out 
the existence of linked configurations at higher Hopf charges. The solutions show both 
similarities and differences to the Nicole and Skyrme-Faddeev models. Solutions clearly 
follow close to a potential lower energy bound of Eafz > Q 3 - Proof of this conjecture 
for both the AFZ and Nicole models would imply a universal feature of Hopf solitons, 
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further to the solutions being knots and (un) links. Technical difficulties required the use 
of a symmetry breaking term to be added to the minimisation which slightly altered the 
problem. This symmetry breaking term was kept small and so the results should give a 
good reflection of the static solutions of the AFZ model. Further numerical or analytical 
evidence would be useful to back up the work undertaken here. 

The symmetry breaking minimisation motivates the investigation of a one parameter 
family of conformal Skyrme-Faddeev models. The two extremal models are the Nicole and 
AFZ models, with linear combinations of the two models completing the set. The family 
of models has been investigated for a Hopf index with qualitatively different topology of 
static solutions in the extremal models. The energy transition as the parameter is varied 
for these models has been found. Details regarding the topological transition of the static 
solutions has been investigated throughout the set of models. There is a smooth change 
in configurations over a short range of 6. It is not clear whether two local minima exist for 
some 6 or not. It is most likely that only one local minima exists for each but that these 
minima become shallow over a small range. The two qualitatively different configurations 
have very similar energies. In fact, the linked configurations look very similar to the axial 
configurations in the transition range. The solution transition is between axial and linked 
configurations. No intermediate configurations exist. The Nicole model contributes more 
energy than the AFZ model for a given coefficient and more than half the models have 
static solutions that are links. In the region where contributions from the Nicole term 
and the AFZ term are roughly equal, static solutions are much harder to find. Volume 
preserving flow finds these minima reasonably well for most models. For parameters near 
to the topological transition of static solutions, solitons are found to be slightly unstable 
up to small perturbations. These solutions appear to have near-zero modes associated 
with the breaking of axial symmetry, at least in the Q = 4 case. A number of issues 
relating to CSF models are left open to investigation. An explicit analysis of the stability 
of solutions to such models would be of use. Do the equations of motion of such models 
simplify in the same way as with the Nicole and AFZ models and can an energy bound 
of the form (26) be proven? 

These results provide further evidence of some universal features of Hopf solitons 
such as the appearance of knotted and linked solutions. Studying models, such as the 
Nicole and AFZ models, can help to explain the structures found in more physical models 
such as the Skyrme-Faddeev model. Volume preserving flow could be used to find static 
solutions to a range scale invariant field theories and might also be useful in the numerical 
construction instantons. 
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